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AIRCRAFT TURBOFAN BOISE 

John F. Groeneweg and Edward J. Rice 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

Recent advances in the understanding of turbofan 
noise generation and suppression in aircraft engines 
are reviewed with particular emphasis on NASA research. 
The review addresses each link in the chain of physical 
processes which connect unsteady flow inte actions with 
fan blades to far field noise. Mechanism identifica- 
tion and description, duct propagation, radiation and 
acoustic suppression are discussed. Recent advances in 
the experimental technique of fan inflow control assure 
that in-flight generation mechanisms are not masked by 
extraneous sources in static tests. Rotor blade sur- 
face pressure and wake velocity measutemenfs aid the 
determination of the types and strengths of the gener- 
ation mechanisms. Approaches to predicting or measur- 
ing acoustic mode content, optimizing treatment impe- 
dance to maximize attenuation, translating impedance 
into porous wall structure and interpreting far field 
directivity patterns are illustrated by comparisons of 
analytical and experimental results. A persistent 
theme of the review is the interdependence of source 
and acoustic treatment design to minimize far field 
noise. Areas requiring further research are discussed 
and the relevance of aircraft turbofan results to 
quieting other turbomachinery installations is 
addressed . 

INTRODUCTION 

Over the past decade the noise generated by the 
fan component of aircraft turbofan engines has been the 
subject of vigorous research. The emphasis on the far. 
reflects the fact that, for the high bypass engines 
which dominate the world fleet of large commercial 
transports, the fan controls flyover noise on landing 
approach and is a strong contributor along with jet 
noise on takeoff. Indeed, for the next generation of 
turbofan engines the prominent contribution of the fan 
to propulsion system noise is projected to continue. 
Figure 1 shows the results of a system noise predic- 
tion done as part of an energy efficient engine design 
study (1). Component noise levels in terms of tone 
corrected perceived noise decibels are given at take- 


off and approach conditions. The fan controls the 
totals at both conditions even with the suppression 
provided by substantial use of acoustic treatment. 
Similar conclusions about the importance of fan noise 
have been drawn in other studies (2,3). The purpose 
of this paper is to review recent results of research 
on fan noise generation and suppression, and to iden- 
tify significant gains and remaining gaps in our under- 
standing and ability to predict and control this annoy- 
ing source. 

The scope of this review is limited to results 
drawn mainly from NASA initiated work carried out in 
roughly the last five years. We believe that this 
definition of scope, relaxed and supplemented in 
specific areas, leads to a reasonably accurate picture 
of the state-of-the-art while not attempting to 
exhaustively cover parallel efforts and results. Sev- 
eral existing reviews provide extensive bibliographies 
and summarize earlier results in turbomachinery noise 
(4), flight effects (5), and suppressors (6,7). Our 
review builds upon, but, primarily, extends and updates 
these earlier efforts to cover significant advances in 
experimental flight simulation techniques, diagnostic 
measurements, theoretical modelling and computation. 

The flow chart in Fig. 2 illustrates the chain 
of physical processes which links unsteady aerodynamics 
of the fan flow field to the resultant far field acous- 
tic signature. Elements in ovals are inputs to, or 
outputs of, the processes in the rectangles. The four 
processes — (i) blade unsteady aerodynamic response, (ii) 
coupling to the duct, (.ii) propagation in the duct 
which may have acoustically treated walls, and (iv) 
acoustic coupling (radiation) to the far field — have 
each been studied and modelled separately as conven- 
ient subdivisions of the overall problem. A knowl- 
edge of the inputs and outputs — (i) unsteady flow 
field disturbances, (ii) blade surface pressure dis- 
i cibutions, (iii) duct acoustic mode content at the 
entrance, and (iv) exit of the duct — is required to 
link the processes and arrive at the final output 
which is far field directivity (and spectra). Of 
course, from an experimental viewpoint, the interme- 
diate inputs or outputs are often missing; only 
acoustic measurements in the far field are available. 
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Figure 1. - Flyover component noise levels for en advanced turt»ofin. 
(Ref. 1) 
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Figure 2. - Elements of fan noise generation and 
suppression. 


In fact, one of the greatest hindrances to applying 
the theories for the Individual processes to practi- 
cal situations has b?en the lack of definition of the 
key Input/output quantities at the Interfaces. Recent 
diagnostic measurements and analyses have begun to 
correct this deficiency. This review Is organized 
around the elements of Fig. 2. 


GENERATION PROCESSES 


Mechanism Identification 

The primary Input to the fan noise generation 
process at eubsonlc tip speeds Is a description of the 
unsteady aerodynamic fluctuations encountered by a 
blade row. In spectral terms, random fluctuations 
produce broadband or narrowband (tone-1 ike) random 
noise while periodic fluctuations produce tones at 


multiples of the blade passing frequency At super- 
sonic tip speeds the added phenomenon of multiple pure 
tones associated with rotor-locked leading edge shocks 
comes Into play at multiples of shoft rotation fre- 
quency. The Initial task of fan ml ie description is 
to Identify the dominant mechanisms in terms of the 
origin of the responsible flow disturbance and the 
blade row with which It Interacts. On the engine 
cross section In Pig. 3. some of the candidate turbo- 
fan mechanisms sre labeled with flow disturbances 
grouped according to the blade row with which they 
Interact. Sample narrowband spectra which Identify 
the components associated with subsonic and super- 
sonic tip speeds sre also shown. Most progress has 
been made In understanding the generation of tones 
which usually dominate the spectrum levels; the rela- 
tive Importance of potential broadband generation 
mechanisms such as those associated with Inlet bound- 
ary layer or wake turbulence remains vague. 

r r." ?! ACOUSTIC TRIATMCN1 



(SUPERSONIC r BlAW PASS,NC 
(SUPERSONIC FWQUtNCY 

r— muitipu 
PURI 
TONIS 



1_L L I 1 

1114 

FRIOUINCV. IHJ 


**TIP 

SUBSONIC 

I 1 — BLAK PASSING riSQUINCY 

r FAN BROADBAND NOISt 



. T..rhnfsn no.ua mnarilinn muhinilllU ifld mCtfl. 


Flow disturbances are divided Into two categories 
In Fig. 2: those originating external to the engine 

but drawn Into the Inlet, and those originating Inside 
the engine. WMle it has long been recognized that 
ingested external disturbances may control tan noise 
generation (8); It was the high bypass engine flyover 
noise data, acquired In connection with noise certifi- 
cation requirements, which established that ground test 
tone levels were controlled by extraneous Inflow dis- 
turbances unrepresentative of flight (5). In fact, 
the practicality of the concept of choosing vane-blade 
ratio for cutoff (9,10) to greatly reduce the funda- 
mental tone was, at first, only confirmed in flight or 
In a wind tunnel as shown by the examples In Fig. 4. 


Flight Simulation . The approach to controlling 
the Inflow for flight simulation In static tests has 
evolved around the concept of Inlet honeycomb-grid 
flow conditioners which must be acoustically trans- 
parent over the frequency range of Interest. Figure 
4 shows the range of Inflow control devices (ICD's) 
Investigated at NASA Lewis (11, 12, i3, 14). The 
sizes of the external devices. Figs. 5(a) and (b) 
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ranged from roughly U to 2 fan diameters. An in- 
duct honeycomb, Fig. 5(c), was found to be aero- 
dynamicalls effective but unacceptable from an 
acoustic transmission standpoint. The first gener- 
ation design, Fig. 5(a), drew on flow conditioning 
work for turbulence reduction (15) to arrive at the 
screen-honeycomb composite structure. The most recent 


version, Fig. 5(b), is reduced in size, uses honeycomb 
only, and employs thinner support ribs with more care- 
fully bonded joints and cleaner attachment to the in- 
let lip. The shape conforms to an equlpotent ial sur- 
face. 

Acoustic results fcr a JT15D engine in terms of 
fundamental blade passing tone and broadband directiv- 
ities are shown in Fig. 6 for the two ICD's. Residual 
tone levels with inflow contvol approach the broadband 
over much of the angle range at these low fan speeds. 
Interestingly, the broadband levels remain essentially 
unchanged with Inflow control indicating that another 
mechanism, probably Internal to the fan, controls this 
spectral component. 

Flight data from the JT15D on an OV-1 test-bed 
aircraft (14) confirm the effectiveness of the ICC of 
Fig. 5(b) as sliown in Fig. 7. The fundamental tone 
directivity with Inflow control agrees well with the 
flight data except at the most forward angles where 
signal-to-noise ratio is low for flight. 

Blade surface pressure measurements have proved 
to be a valuable diagnostic tool in evaluating the 
quality of inflow tc the fan and, with Inflow control, 
in determining the residual internal sources control- 
ling flight levels. Miniature pressure transducers 
mounted near the fan blade leading and trailing edges 
at several spanwlse locations are used in conjunction 
with a telemetry system (8,12,16) to continuously sur- 
vey the circumferential variation of unsteady blade 
pressures. This technique originally identified long- 
itudinally persistent, circumferentially localized 
disturbances attributed to atmospheric turbulence 
elongated by the streamtube contraction in the inflow 
(8). Such disturbances, which may also be caused by 
Ingested vortices, wakes and instabilities associated 
with flow around the inlet lip, produce scrong narrow- 
band random tones. 
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Figure 5. - Inflow control devices Air 'light fin noise simulition tested it Lewis Research Center. iRef. 121 
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Figure 6. - Far field directivity with inflow control devices in JT15D 
engine. Af ■ 25 Hz. (Ref. 12) 



Figure 7. - Effectiveness of inflow control for flight fan tone 
simulation. (Ref. 141 


Figure 8 compares narrowband blade pressure spec- 
tra without and with Inflow control. Without an ICD, 
Fig. 8(a), the spectrum shows strong harmonic content 
at all multiples of shaft rotation frequency resulting 
from multiple encounters of the blade transducer with 
circumferentially varying flow disturbances. The 
additional scales on the abscissa are distortion number 
(multiple of shaft frequency) and the circumferential 
acoust'.c mode number corresponding to blade number 
minus distortion number. Inflow control eliminates the 
randomly varying disturbances and the corresponding 
bulk of the shaft harmonics as illustrated in rig. 

8(b). Those distortion numbers that remain are assoc- 
iated with periodic, internally generated flow distur- 
bances which are fixed in space or have fixed rotation 
rates with respect to the rotor. As a result, clues 
to the mechanisms governing flight levels are found 
from the prominent residual peaks. 
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Figure 9. - Directivity pattern indicating presence of 22-lobed circum- 
ferential mode due to rotor-strut interaction, JT15D fan, 10 500 rpm. 
(Ref 12.1 


Substantial effort has also been applied to the 
inflow control problem by industry (16-23) Including 
flycver level comparisons to static projections (22) 
and development of ICD design procedures (23) . The 
first generation of large engine ICD's, roughly 3 fan 
diameters in size, is currently in use. Although the 
quantitative agreement of inflow control and flight 
data is still subject to some improvement, the current 
"state-of-the-art" does allow the study of bona fide 
internal sources controlling fan noise generation in 
flight. An alternative to ICD's are anechoic wind 
tunnels ( 24-’’6) which also have been found to eliminate 
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the bulk of the extraneous Inlet disturbances 

In-Flight Sources . Once the study of Internal 
mechanisms Is made possible by Inflow control, the 
task becomes one of Identifying the Interactions re- 
sponsible for the tone levels observed over the range 
of engine speeds. Rotor wake-Btator interaction re- 
mains a prime mechanism; hut, *’>*n with thn bludu -vaim 
ratio chosen to prevent fundamental tone propagation, 
other Interactions may come Into play. For example, 
the JT15D engine exhibits a strong fan fundamental 
tone which appears at a speed corresponding to the 
start of propagation of the 22-lobed acoustic mode as 
shown in Fig. 9(12). The source of the 22-lobed 
acoustic mode Is the Interaction between the 28 fan 
blades and the six structural support struts down- 
stream of the fan stator. The blade pressure spectrum 
in Fig. 8(b) shows that a strong 6-per revolution dis- 
turbance Is sensed on the rotor. A prime candidate 
for the interaction mechanism is a strut potential 
field extending upstream through the stators and 
interacting with the rotor. An alternate explanation 
would be the interaction of residual rotor wakes with 
the six engine struts generating the 22-lobed spinning 
acoustic mode which is sensed on the rotor as a 6-per 
revolution disturbance. Existing large high-bypass 
turbofans also contain downstream struts. The next 
generation of engines will incorporate integral strut- 
stator vane assemblies with a potential for still more 
complicated interactions (27). 


^ENGINE (STDIMNE 



Figure R - Variation of broadband levels with (an operating 
point, JT150 (an in anechoic chamber Af • 80 Hz. 

10500 rpm. 

As mentioned earlier, the random disturbances re- 
sponsible for broadband generation appear to be inter- 
nal. Broadband levels vary strongly with fan operat- 
ing point (rotor incidence angle or loading) as shown 
in Fig. 10. Fan blade suction surface flow separation 
and interaction with the trailing edge, blade tip 
interaction with the casing boundary layer and rotor 
wake turbulence (midspan or tip) interaction with the 
stator are candidate mechanisms although the latter 
seems to be discounted by rotor alone experiments (28). 

While the multiple pure tone (MPT) generation 
mechanism is clearly associated with the rotor leading 
edge shocks and their blade-to-blade nonuniformity, 
quantitative descriptions of the source are lacking. 

Two aspects of the process have been investigated in 
connection with concepts to reduce this strong inlet 
source at takeoff fan speeds. First, a radical fan 
with compound leading edge sweep was designed to keep 
the normal component of blade inlet relative Hach num- 


ber subsonic over the entire span (30). Except for 
blade end effects and the sweep reversal point, a 
major portion o f the strong leading edge shock system 
was expected to be eliminated. Figure 11 shows the 
measured MPT power results obtained with the swept 
design compareJ to a conventional, unswept fan (31). 
Sweep delayed .he onset of MPT's and reduced the levels 
over s portion of the tip-speed range including 

takeoff. The second aspect of inlet noise generation 
at supersonic tip speeds concerned the observation 
that total tone power peaks beyond the transonic speed 
and then falls off. A fan designed for unusually high 
specific flow, 220kp/sec .m^ , at high tip speed, 553 m/ 
sec, exhibited a marked tone power decrease at design 
(32) although not qualitatively different from other 
high tip speed designs. The phenomenon appears to be 
only partially attributable to propagation inhibiting 
effects of elevated inlet Mach numbers and may be 
associated with non-linear propagation characteristics 
in combination with the angle and associated strength 
variations of the leading edge shocks (33). 



HP INLET RELATIVE MACH NO. 

Figure 11. - Effect o( rotor leading edge sweep on multiple pure tone 
generation. iRe(. ill 


Mechanism Desc r iption 


Ducted C asca de Response . Considerable effort has 
been expended within the last five years to model the 
noncompact compressible response of a ducted cascade 
of blades to unsteady upwash velocities. Perhaps the 
most complete description available is the three- 
dimensional lifting surface theory (34) for a rotating 
cascade in an annular duct. This blade response and 
duct coupling analysis is the heart of specialized 
studies of rotor-inflow distortion (35) and rotor- 
stator interaction (36). These linear analyses are 
for the dipole-type sources at the surface of a cas- 
case of thin (in some cases twisted) blades and repre- 
sent exact solutions to the linearllzed continuity and 
momentum equations (37, Chap. 5). 

Three features of these analyses are considered to 
be important advances. First, the 3D approach makes 
possible the calculation of the circumferential and 
radial content of acoustic modes in annular or cylin- 
drical ducts; the complete description of modal content 
is precisely the input required for successive propa- 
gation analyses. Second, cascade analysis predicts 
chordwise unsteady pressure distributions and inte- 
grated responses which differ substantially from sin- 
gle blade results (38,39) that ignore blade-to-blade 
interactions (solidity) and the interblade phase angle 
of the disturbance. Third, source non-compactness, 
retained by calculating chordwise in addition to span- 
wise pressure fluctuations, has been shown to produce 
significant differences in calculated power compared 
to compact analyses. The magnitudes of the differ- 
ences, which depend on incident disturbance shape and 


propagation direction with raipcct to the scan flow, 
are hlgheet for elngle distortion modes (38). For 
realistic distortion profiles represented by a conbln- 
atlon of distortion nodes, the effects of nonconpact- 
ness are less dramatic with a tendency for the compact 
analysis to underestimate fundamental tone power for 
upstream propagation and overestimate the power propa- 
gating downs tream( 35) . 

Since full 3D calculations are complex and lengthy, 
• quael-3D analysis, which uses 2-D (strip) theory tor 
aerodynamic responae but annular duct acoustics for 
modal prediction, was Investigated (35). The results 
Indicated that the quasl-3D approach produced rela- 
tively small errors In power, greatly reduced compu- 
tation time, and fulfilled the requirement to predict 
annular duct acoustic modes. Consequently, the quesl- 
3D approach was adopted In the development of a compu- 
ter program (40) which considered three types of flow 
disturbances: Inlet turbulence, rotor mean wakes and 

rotor wake turbulence. This quasl-3D approach Is still 
being evaluated by data-theory comparisons. 

The 3-D lifting surface tone power predictions 
have been compared to fan noise data (41) where the 
controlled source consisted of the fan Interacting with 
an array of Inlet distortion rod wakes. Figure 12 
shows excellent agreement between the predicted total 
Inlet fundamental tone power as a function of far. speed 
and the measured narrowband tone power obtained from 
far-fleld measurements. Note the changing mix of 
radial mode contributions to the totals and the non- 
monotonic increase with speed in both theory and data. 



6500 7500 8500 9500 10500 11500 


FAN SPIED, rpm 

fi 9 ur* 12. - Speed fcpenfcnce of inlet fundementil lone 
modal power qeneraled by 41 rod wakes interacting 
with a JT1 50 fan. iRef. 41i 

The Intermediate quantity between blade response 
and duct coupling is blade pressure (Fig. 2). The 
cascade response portion of the code in (40), was used 
to calculate the chordwlse magnitude of the unsteady 
blade pressures due to Interaction with Gaussian wakes 
produced by upstream radial rods. As shown in Fig. 13, 
the high disturbance frequency associated with many 
s 4 1 ) rod wakes is predicted to produce many rapid 
changes In pressure along the chord. Typical miniature 
transducer sizes are Indicated near the leading and 
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trailing edges. For high disturbance frequencies the 
analysis Indicates that measured blade pressure ampli- 
tudes are subject to uncertainty due to finite trans- 
ducer sire and sensitivity to transducer location. 
However, experimental checks of the cascade response 
analysis using carefully controlled flow disturbances 
are needed. 

Rotor - Stator Interaction . To apply the cascade 
response analyses just described to one of the main 
tone noise generation mechanisms, rotor-stator inter- 
action, a thorough description of the rotor produced 
disturbance flow field Is required. The need to des- 
cribe blade wakes has long been recognized and a large 
body of wake data Including mean and turbulence proper- 
ties has been accumulated on laboratory fans (e.g.. 42, 
43). In addition to mid-span wakes, secondary flows 
such as tip vortices have been recognized as potential 
noise contributors (44). Therefore, a linear cascade 
analysis Including spanwise gust components has been 
developed to allow the relative noise contributionr of 
tip vortices and mid-span wakes to be determined (45). 
What is lacking is a thorough model of the total rotor 
downstream flow field which Is linked to fan design 
parameters and Is validated by experimental data. 

Some wake data have been obtained as functions of 
downstream distance on a fan operated with forward 
velocity In an anecholc wind tunnel (46). Mean wake 
stator upwash velocity profiles are shown in Fig. 14 
as a function of spanwise position. The magnitudes 
vary substantially with radial location, but most sig- 
nificantly the profile near the tip Is characterized 
by an extra upwash cycle between successive blades 
corresponding to strong secondary flows, probably a 
tip vortex. The variation of stator upwash harmonics, 
the required Input to generation analyses, are shown 
In Fig. 15 as a function of downstream distance. From 
the complex variations observed, it must be concluded 
that simple Gaussian profiles which decay and spread 
monotonlcally with distance are an Inadequate des- 
cription of this flow field. 

Acoustic data are available from rotor-stator 
spacing experiments on the same fan as was used for 
the wake measurements just described. Two stator vane- 
rotor blade ratios were examined: one for propagating 

and the other for cutoff fundamental tones. Figure 16 
shows the Inlet narrowband tone harmonic power level 
variation with rotor-stator spacing. Residual ’cvels 
of the fundamental for the cutoff case '25 vanes) are 
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Figure I A * Rotor mean wake velocity profiles as a func- 
tion of spanwise location. 1. 23 rotor chords downstream. 
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DOWNSTREAM DISTANCE tROTOR BLADE CHORDS! 
tai 9% of span from tip. 

(b) 30H of span from tip. 

Fiqure 15. - Mean wake harmonic levels from ensem- 
ble average spectra. 80* design speed, upwash com- 
ponent. tunnel vel. ■ 41 m/sec. tRef. 46) 


80* design speed, tunnel velocity ■ 41 m/sec. (Ret. 46) 

nearly constant suggesting that a weak Inflow distur- 
bance-rotor interaction governs in this case rather 
than a rotor-stator interaction due to stator blade 
nor. -uniformities sufficient to generate other propa- 
gating modes (48). Note that 25-vane second and third 
harmonic levels are higher than corresponding harmonic 
generated by the 11-vane set Indicating a difference 
in the response and/or coupling to acoustic modes of 
the two stetors. The 11-vane stater had longer chords 
than the 25-vane stator in order to maintain the same 
solidity. Comparisons between the measurements and 
theory in Ref. 40 are in progress using the corres- 



ponding wake measurements as input. Tone powers meas- 


ured in rotor-stator spacing and vane-blade ratio 
experiments in an anecholc chamber using inflowcontrol 
have been compared to a 2-D (strip) model with encour- 
aging results (49). Wake data were not acquired so a 
wake model was used. While 2-D theory may do rela- 
tively well for power predictions; calculating far 
field directivity and, therefore, acoustic mode con- 
tent requires more sophistication in handling duct 
geometry and, probably, in describing the wake/vortex 
flow field. 

Source Modal Content . Both the magnitudes and 
phases of all the acoustic modes generated are the 
fundamental Inputs to propagation analyses. To date, 
most available mode information is much less complete. 
For tone sources circumferential mode numbers may be 
determined from rules about blade row interactions 
(50). No such simple rules exist for radial mode 
numbers. For example, in rotor-stator Interactions the 
rotor wakes become increasingly radially skewed with 
downstream distance; one wake may simultaneously 



Figure 16. • Viriition of Mrrowbind tone h»i ironic power 
levels with rotor-stltor spicing, 80* design speed, 15 
rotor blides. tunnel velocity • 41 m'sec. (Ref. 47. ' 
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intersect several blades with the intersection points 
sweeping radially with time. In multimodal situations 
corresponding to random flow excitation, all possible 
modes that the duct can support are candidates. 

Experimentally, two approaches to mode identifi- 
cation have been used. In the first, modal content 
is Inferred from far field directivity. Individual 
principal mode ahapes obtained from radiation theory 
are matched to the observed lobe shapes (51,52). This 
method ignores any phasing effects between modes in 
the far field, depends on the accuracy of the radia- 
tion model and breaks down for lnolvldual modes as the 
number of propagating modes Increases. In the limit 
where a large number of modes propagate and a smooth 
directivity pattern results, either equal amplitude or 
equal energy modal distributions have been assumed and 
compared to data. One auch study indicated that the 
equal energy assumption produced a better fit to inlet 
tone directivities than equal amplitude (S3); in 
another study (54) the equal energy assumption gave 
fair fits to the multimodal directivities in some 
cases but poorer results in others, particularly at 
aft angles. 

The second experlL. .tel approach has Involved 
attempts at direct mode determination from ln-duct 
microphone measurements. Techniques range from matrix 
inversion of N selected wall pressures measurements to 
deterministically solve for N preselected modes (55), 
to a least-squares approach where the number of mea- 
surements is at least twice the number of modes (56). 
Formidable practical difficulties exist: radial mea- 
suremznts upstream of the fan Introduce distortions 
and their associated extraneous modes, and measure- 
ments on the wall alone require large numbers of 
microphones distributed axially and circumferentially. 
Demonstrations of the techniques have been United to 
relatively low speed fans with conditions rather far 
removed from the turbofans of Interest. Direct 
measurement techniques require additional development 
and still fall in the category of research efforts, 
not routine tools. Thus, results from 3-D or quasi 
3-D analyses (35, 36, 40, 41) are presently the main 
source of detailed modal information. 

PROPAGATION AND RADIATION - UNLINED DUCTS 

In this section we de«l with propagation in and 
radiation from ducts with hard walls. Untreated cases 
are important for carrying generation predictions to 
the far field to compare with experiments. Treated 
ducts will be discussed in the next section. 

Blade Row Transmission 

Sound generated by pressure fluctuations on the 
stator can only propagate upstream by traversing the 
rotor and vice-versa. Inlet levels of propagating 
rotor-stator interaction tones sometimes display low 
sensitivity to rotor-stator spacing (51, 57) indicat- 
ing substantial rotor attenuation. Stator generated 
modes incident on the rotor may be acattered into 
modes at other blade passing harmonics (58) or sum and 
difference tones in multistage arrangements (59), and 
are more highly attenuated if they are rotating oppo- 
site to the rotor. While actuator disk (60) and 
cascade analyses (61) exist, quantitative experimental 
verification has been limited to two-dimensional 
approximations (62). The existing analytical frame- 
work remains to be applied to 3-D transmission in 
annular ducts for practical cases where the tone wave- 
lengths of interest are often equal to or less than 
the blade chords. 

Inlet Radiation 

Assuming any blade row transmission effects have 
been accounted for, the modal content propagates to 


the far field through variable area ducts carrying 
flow and radiates from duct openings through non-uni- 
form flow fields. For inlet radiation two distinct 
flow fields are of Interest: largely radial potent'al 

flow in the static test case, and an inlet stream tube 
only slightly larger than the inlet diameter in the 
flight case. The Wlener-Hopf technique, ippllcable 
only to inlet lips of negligible thickness, has been 
applied to two idealized cases. One is uniform exter- 
nal and internal flow at the same Mach number, and lh. 
other is a constant Mach number external flow bounding 
a cylinder of higher uniform Mach number extending out 
of the inlet. The former is an approximation to the 
flight case but the latter is unrepresentative of any 
real inlet flow. Two other approaches to analyzing 
inlet radiation have been followed. The first uses 
simplifying assumptions based on ray acoustic while th< 
second uses a fully numerical solution incorporating 
the actual flow field and inlet lip geometry. 

Ray Acoustics in Terms of Cutoff Ratio . Approx- 
imate expressions for inlet radiation have been devel - 
oped in terms of mode cutoff ratio £ * nri/ a ra n \ 1 -‘■f,-’" 
where the dimensionless frequency n “ f D/c; with the 
eigenvalue of the (m,n) tf> mode.amn: duct diameter, U; 
duct Mach number, Mr,; and epeed of sound, c. The key 
simplification realized by the cutoff ratio formula- 
tion is that modes with the same £ and r) propagate 
similarly to the far field. This has been demonstrat- 
ed for radiation from a flanged duct without flow (63) 
and is fairly accurate for principal lobe radiation 
(64). Two important duct mode propagation angles, 
and ii< x , are defined in (65) as 


cos 


and 

cos 1 * 1 , 


where 

S - ^ 



1/(2 


(1) 


( 2 ) 


(3) 


Here, $ and i|/ are the angles which the vector no-mal 
to the wave front and the group velocity vector, re- 
spectively, make with the duct axis. The duct mode 
angle ik, , given by Eq . (2), has been shown to closely 
approximate the angular location of the principal lobe 
in the tar field (65). This conclusion was reached by- 
inspection of the directivity coefficient appearing in 
the Wiener-Hopf solution for the case of uniform flow 
everywhere (66); an expression for the principal lobe 
angle identical to Eq . (£) was obtained. The approx- 
imate equality of ln-duct propagation angle and far 
field principal lobe radiation angle suggests that 
ray acoustics arguments can be used to link the two 
angles for cases where the flow is not uniform. 

For example, ray acoustics ideas have been 
applied to the case where far field velocity is sub- 
stantially less than inlet duct velocity, the limit 
being the static case where far field velocity is 
zero. Based on a ray acoustics analysis which showed 
that refraction in a potential flow is of second order 
in Mach number (67), the wavefronts were assumed to be 
unbent going from duct to far field. That is, was 
assumed to be unchanged. Since $ and 1 are identi- 
cal if Mach number is zero, the giwup velocity in the 
far field was assumed to have been shifted. At M - 
-0.4 and f, % 1 (near cutoff), the calculated radia- 
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tlon peak la at 66° while the group velocity in the 
duct propagates at * 90 . A peak near 66 3 was ob- 
served in the far field fora nearly cutoff aode gener- 
ated by a controlled fan source (52). However, the 
agreement of this observed peak with the theory which 
neglects lip shape may be misleading. A propagation 
phenomenon associated with the very thick inlet lip 
used in the experiment may have controlled the princi- 
pal lobe location. An analysis of propagation in a 
variable area duct with gentle area variation shows 
that mode identity is preserved (68, 69) (i.e., no 
scattering occurs). Thus, as a mode propagates from 
the ir.let throat to the highlight, C increases caus- 
ing and i)' a to decrease. Recent extensions of ray 
theory for propagation through an lrrotatlonal flow 
(70, 71) imply that it is the group velocity vector 
which is unchanged, not the normal to the wavefronts. 
The difference between the two assumptions is signifi- 
cant; e.g.,66 vs 90 peak near cutoff, and current 
evidence points to preservation of group velocity as 
the better approach. Controlled experiments and possi- 
bly numerical simulations are needed to settle this 
issue . 

Numerical Model . A hybrid numerical program has 
been developed (72) and exercised (73) to calculate 
both the internal and external sound propagation for 
actual engine inlet geometry ar.d flow conditions. It 
Is a hybrid program in the sense that a finite element 
method is used within the duct and in the near field 
and an integral radiation method handles the far field. 
Iteration is tequlred to match the two solutions at the 
interface. A potential flow program is used to gener- 
ate the steady flow for the actual inlet geometry; 
boundary layers are not included. The input to the 
program is the pressure profile for a given mode in the 
annulus at the fan source. While the program cannot 
handle the combination of high Mach number and high 
frequency due to computer storage limitations, some 
inlet geometry effects have been studied which were 
previously impossible to analyte. 

Figure 17 compares the numerically predicted inlet 
tone directivity to the measured levels generated by 
the same controlled source mentioned previously — a 
JT15D engine with inlet rods (52). A single (13,0) 
mode propagates at the fan speed shown. The excellent 
agreement between the hybrid solution and the data is 
in contrast to the Wlener-Hopf solution for an infin- 
itely thin lip, which is also shown. The thick lip 
used in the experiment ( thlckness-t o -d iamet er ratio of 
0.5) shifts the radiation peek toward the axis, as 
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Figure IT. - Sinqle mode inlet directivity companion ot 
theory with e«per:ment 31X1 Ht. il 3, O' mode. iR*f. 86i 


discussed in tbs proceeding section, snd sets ss a 
shield to reduce the levels in the aft quadrant. The 
dependence of the directivity on inlet llo thickness is 
illustrated in Fig. 18 where the shielding effect is 
also clearly evident. The numerical results show that 
the radiation peak moves aft ss the lip gets thinner. 

At s thlckness-to-dlameter ratio of 0.1, the rsdiatlon 
pattern agrees very well with the Wlener-Hopf (zero 
thickness) result shown in Fig. 17. The hybrid pro- 
gram is a powerful tool for the solution of "real- 
world" inlet rsdiatlon problems. 



ANGIE FROM INLET, deg 

figure 18. - Eltect al inlet lip thick- 
ness on single mode directivity. 
M D ■ a 3150 Hz (13.01 mode, 
(theory normalized to 100 dB at 
HP l (Ret. 85) 


Exhaust Rad iatio n 

In contrast to the complex inlet flow field, the 
exhaust flow, neglecting mixing, is much simpler. The 
fan exhaust may be approximated as an emerging cylin- 
drical flow at Mach number M^ surrounded by a uniform 
flow at Mach number M^ which fits the requirements for 
an exact Wlener-Hopf radiation solution. The ray 
acoustics, mode cutoff r„t ir approach to an approxi- 
mate solution can also be applied with more confidence 
to the aft slip layer. Starting from the zero-flow 
flanged duct radiation equation, a coordinate trans- 
formation was applied to account for the duct flow, 
and ray acoustics arguments were applied across the 
slip laver (74). Single mode aft directivity from the 
approximate expression is compared to the full Wlener- 
Hopf solution (75) in Fig. 19. The good agreement 
builds confidence in the simplifications used to gen- 
erate the approximate solution. The Wiener-Hopf solu- 
tion gives finite levels in the zone of silence in 
contrast to the ray acoustics result although the par- 
ticular values from (75) are believed to be incorrect. 

The location of the principal lobe in the far 
field. ,, is found from the approximate theory (74) 
to be: ,p _ 

cos , fp . - V f ->P-l (1 -Mp 2 ) (4) 

<‘- m d 2) <*• +m dv ,C2 - 1) 


for the static case (M • 0). For • 4-0 f< and £ ■ 
I, v. ■ 160° measured from the exhaust axis indicat- 
ing iRjt modes near cutoff radiate to the inlet quad- 
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rant. The analogous Inlet analysis , Eq. (2), 
indicated that near cutoff node peaks remain In the 
Inlet quadrant. Thus, the Inlet qusdrant contslns the 
near cutoff node peaks no matter where the sound 
originated. 


SUPPRESSION - LINED DUCTS 


Ducts between the turbomachinery components and 
the observer can be lined with sound absorbing mater- 
ials to greatly reduce the radiated noise. Early 
workers (77) on suppressors considered the use of 
splitter rings In the duct to Increase treated area 
and decrease the distance between treated surfaces. 

Current emphasis Is on the optimization of wall 
treatment alone to minimize aerodynamic losses and 
weight penalties In aircraft applications. The atten- 
uation achieved is very sensitive to the source modal 
characteristics used as Input (78, 79). Input cases 
of interest range from a limited number of modes asso- 
ciated with tone noise from periodic blsde row Inter- 
actions, to multi-modal situations associated with 
random procasses exciting all the modes the duct can 
support. The latter number may be very large. At 
high frequencies considering modes spinning In one 
direction, the number can be estimated from ref. (80) as 
N ■ T 2 n 2 /8(1 - M_ 2 ) . For typical values of » -0.4 
and n * 20 at blade passing frequency, N*£>00. In 
such esses, some method of handling the modal distri- 
bution as a continuum subject to a simple rule des- 
cribing the energy distribution Is desirable. 

The status of three aspects of suppressor re- 
search will be discussed: analytical propagation 

approaches Including the cutoff rat.j method, numer- 
ical propagation programs, and suppressor materials 
characterisation . 


Duct Acoustics - Analytics'. 

Much analytical effort has gone Into describing 
propagation of acouitlc modes In simple geometries 
such as cylindrical and rectangular ducts lined with 
sound absorbing eaterlals. A complex acoustic Impe- 
dance (reslstan e and reactance) Is used to specify 
the wall boundary condition for point reacting liners. 
Solution of l ne wave equation in the duct results in a 
complex eigenvalue and wave number for each mode. The 
real part o. the wave number defines the attenuation 
while the imaginary part defines the axial wave speed. 

From such single mode solutions, attenustlon con- 
tours plotted In the lmped-nce plane, as shown In Fig. 
20, reveal an optimum Impedance for which the maximum 
possible attenustlon Is obtained. The Imperil yt para- 
meters are duct Mach masher, boundary layer thickness 
and frequency. Srudles of optimum impedance led to 
the discovery (81) that mode cutoff ratio was a key 


♦ MAXIMUM ATTENUATION DEFINING 
OPTIMUM IMPEDANCE 



Figure 20. - Sound ittenuitior contours boundiry layer 
thickness t • 0; liner length/passage height, UH • 1.0. 
IRef. 82.1 



Fiqure fl. - Opt'mcm mall rasistance as a junction o! 
mode cutoff ratio lor various inlet wall h mndsry 
layer thicknesses. iRef. 8?i 

correlating parameter: modes with similar cutoff 

ratios propagate In a similar fashion. Figure 21 
shows the correlation of optimum resistance as a 
function of cutoff ratio for different boundary layer 
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thicknesses. Each point rcpraaanta a tingle node 
calculation. The curvet through the point* are from a 
correlation equation reported In Ref. (82) which 
alao contalna correlatlone for maximum poaalble atten- 
uation and optimum reactance. Hlgheet attenuation 
occurs near cutoff and drops rapidly with Increasing 
cutoff ratio. Cutoff ratio correlations have bean 
found for far field directivity (63) duct termination 
loss (63) and transmission loss through a variable 
area duct based on t at the throat (68, 69). 

Multimodal casea, such as the example of N * 600 
cited above, ran be handled as a continuum in cutoff 
ratio. The modal masher density function In a duct 
la Riven In (80) by: 


1 dN 2 
N dC * 7? 


(5) 


The number density Is converted Into a modal power 
density by multiplying by a weighting function which 
must then be estimated by a technique such as a least 
squares fit to experimental hardwall directivity pat- 
terns as described In Ref. (6). 

The goal of these approximate suppressor analyses 
is to predict the suppressed far field directivity. 
Additional work Is needed to Improve the quantitative 
results by Including modal scattering at the hard-soft 
ll'ier Interface, refraction around the inlet lip at 
hi;h Mach number and redirection of sound by the inlet 
lip contour. Analogous optimum Impedance and attenu- 
ation correlations remain to be derived for the aft 
duct case. 


Duct Acoustics - Numerical 

Recent reviews oi the application of numerical 
methods to iuct acoustics appear In Refs. (83) and (84). 
The hybrid numerical program in the "Inlet Radiation" 
section also handles the case of acoustically treated 
Inlet walls. Comparisons of the calculated and meas- 
ured suppression for a series of very short Inlet 
liners (L/D»0.15) are shown In Fig. 22 (85). The 
experiment again used a JT15D engine configured to 



ftquri 72 . - Companion of numeuwl suppression/ 
radiation analyst! to experiment as a function of 
liner resistance il 3. 0> mode, 3150 Ht, * ■ 0. 5. iR»f. 
85. 


produce a single (13,0) mode. Three different liner 
resistances were tested (86). For a single mode the 
attenuation la independent of angle. Therefore, at 
each resistance, data points are plotted for angles 
between 50° and 80° where the single mode at blade 
passing frequency dominates the levels. The calculated 
values agree well with the data which vary less than 
5 dB with angle. 

The hybrid program has been used to calculate 
attenuation contours in the Impedance plane and thus 
define the optimum Impedance (85) for the single mode 
as shown In Fig. 23. Optimum values of resistance 
and reactance are 0.6 and 0.85 compared to values of 
1.14 and 0.5 calculated for a single (13,0) mode in 
the soft-walled duct. Corresponding ma lmum possible 
attenuations differ by 20 dB (50 dB numerical vs. 30 
dB single mode). These comparisons snow the Impor- 
tance of the Inclusion In the hybrid program of modal 
si itterlr.g at the hard-soft interface, particularly 
for a near cutoff mode entering a short liner. As the 
liner is lengthened, the numerical calculation 
approaches the analytical result t-hich neglects scat- 
tering. A semi-emplrlcal correction was applied In 
Ref. (85), but explicit Inclusion of scattering in the 
analysis is preferable. 

The challenge In applying the numerical approach 
is to remove the limitations on frequency and Mach 
number ranges Imposed by computer storage and run 
time. One possibility is to develop a transient sol- 
ution method which can potentially cut storage re- 
quirements by several orders of magnitude (87, 88). 
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Fiqure 23, - Calculated optimum impedances from a hybrid numerical 
solution with scattering compared to the single soft «ali mode solu- 
tion. 3)50 Ht (13.0) mode, liner l/dj, • .15. iRef. !5i 

Suppressor Mat erials 

In order to realize the attenuations calculated 
bv the propagation analyses, the duct boundary condi- 
tions expressed as acoustic Impedance must be trans- 
lated Into liner construction parameters, e.g-, treat- 
ment thickness, wall porosity, fihe' sizes and bulk 
densities, etc. For extended reaction liners, normal 
wall Impedance is Inappropriate and coefficients in 
the wave equation describing propagation in the liner 
Itself must be related to real bulk absorbing material 
properties. The following discussion will briefly 
summarize the physics of the various dissipation mech- 
anisms involved and some recent modelling efforts to 
link the suppressor construction details to the global 
properties needed for propagation analysis. 
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Impedance - Poin t Reecting Materiel *. Helmholtz 
resonator array* formed by bonding a perforated plate 
to honeycomb have been widely ueed to acouatlcally 
treat aircraft engine duct*. In the abaenc* of the 
grazing flow preeent In engine application*, these 
suppressor* are extremely non-linear, i.e. re.itienc* 
depends on sound pressure level because the oscilla- 
tory orifice flow Is govern'd by the Bernoulli equa- 
tion. However, at high grazing flows the resistance lr 
prop it clonal to free stream M/aeh number and Is Inde- 
pendent of sound pressure level. The control of 
resistance by grezlng flow Is due to orifice blockage 
which occurs as the mean flow Is periodically turned 
Into the orifice* and >.h*n separates, reducing the 
discharge coefficients or the effective open area 
ratio of the perforated plate (89). Acoustic diss. - 
pat Ion occurs as the kinetic energy of the hydrody- 
namic Jet-like flow Induced In the o. if Ices Is dissi- 
pated In the surrounding fluid with little pressure 
recovery. Actually, boundary layer properties must 
also be considered sine, the orifice* Ingest flow from 
distances of the order ol one orfflce diameter. 
Acoustic reslslsnce models have been published (66,90) 
as s function of free stream Mach number sad boundary 
layer thickness based on s l/7th powet profile. Re- 
cently, a model using friction velocity as the con- 
trolling parameter has been formulated (91) to handle 
the more general case of a boundary layer ’hlckened 
by diffusion or altered In a way which divorces near 
wall velocities from free stream conditions. 

Perforated p’ate-honeycomb treatment has cost and 
load carrying advantages , but Its acoustic design Is 
complicated by the requirement to know local flow con- 
ditions which vary with engine speed. (The associated 
resistance variation Is favorable In an inlet but det- 
rimental In the exhaust). To overcome the sensitivity 
of resistance to flow conditions and maintain acoustic 
linearity, treated surfaces having very small pores 
formed by densely packed wires or extremely fine 
sereins have been used. In this case, the dissipation 
mechanism is viscous pressure loss In .he fine pores 
as nearly stagnant fluid very near the wall is Inges- 
ted. An impedance model for square weave screens sup- 
ported by perforated pi: •* has been developed (92). 
Mechanical problems associated with fine screens have 
beer ''dressed. Corrosion at the bonds wltn dissim- 
ilar metals can be prevented with coatings. Contam- 
lmation is prevented by the sel f-cleaning action of 
the oscillatory velocities induced in the pores. Aside 
from the increased in situ Impedance predictability, 
fine-pored materials have no acoustic advantage over 
perforated plate when used In a honeycomb-bad rd wall 
treatment . 

Bul k Absorbers . Extended reaction liners formed 
by an uninterrupted layer of densely packed fibers 
retained bv a high porosity face sheet do offer def- 
inite acoustic advantages. The distributed dissipa- 
tion damps back cavity resonances which smooths the 
reactance as a function of frequency and increases 
high frequency absorption. 

At low frequencies the speed of sound in the 
fiber jus matrix is decreasec giving a greater apparent 
treatment depth and enhancing low frequency absorption. 
Beyond these inherently point reading arguments, the 
possibility that axial and circumferential wave travel 
In the material may be used to advantage is being in- 
vestigated (93). Studies of the effects of high sound 
intensities on bulk absorber properties have beer, com- 
pleted (94) and aodels for a range of bulk absorber 
types are under development (95). 

There are mechanical difficulties associated with 


extend'd reaction llnern which must be overcome. Ab- 
sorption of liquids will destroy the acoustic effec- 
tiveness and can pose a safety problem. Non-wetting 
coatings may solve the problem for some liquids and 
allow removal of others using compatible solvents. To 
provide structural strength while preserving extended 
reectlon, porous (acoustically-transmitting) axial and 
circumferential face plate supports may be required. 

CONCLUDING REMARKS 

Status and Outlook 

This review has addressed the links in the -haln 
of aeroacouxtlc processes which c meet turbofan noise 
generation and suppression to the observed far field 
noise. Although the theoretical outlines and a body 
of often flawed empirical data have been In hand for a 
decade or more, the unified application of theory and 
controlled experiment to practical cases Is just now 
occurring. This productive development was delayed by- 
contaminated data and theoretical complexity whose 
physical imp 1 lca r ions were Inadequately communicated 
to or assimilated by applied noise researchers. During 
the part half uecede, key advances nave been made In 
experimental technique and theoretical application 
which open the way to a much b" -lader understanding and 
conttcl of turbofan noise. 

With respect to generation, the development of 
effective inflow control techniques makes possible the 
conduct of definitive experiments on Internally 
controlled blade row interactions. The Initial round 
of such experiments, focused on rotor-stator tone gen- 
eration, has highlighted the question of Lne relative 
importance of secondary flow disturbances (c.g., tip 
vortices) compared to mid-span wakes. Both Improved 
descriptions of the rotor-produced flow disturbances 
and the noise generation computer codes to use them 
ire being developed to answer this question. With the 
help of rotor blade pressure measurements at.d Inflow 
control, rotor interaction with struts downstream oi 
the stator has been identified as a significant noise 
source. The experimental tools now exist to uncover 
ntf.-r such mechanisms which set tone levels In speci- 
fic engines. Computer codes which Incorporate non- 
compact, cascade response to calculate the generated 
acoustic mode content are in th» process of being val- 
idated by experiment. Tone power comparisons show 
good agreement. The next level of validation involves 
far field directivity which depends on a prediction of 
individual modes and their propagation behavior. At 
least a quasl-3D calculation approach must be used to 
give both circumferential and radial m.<dt content. 

Perdlng the development of well-proven modal 
measurement techniques, confirmation cf predicted mode 
content is Intertwined with hardwall propagation and 
radiation analysis. The throret teal approaches dis- 
cussed include approximate analyses. Wlener-Hopf lech 
nlques and direct numerical solutions. Geometric 
acoustics approximations cast in terms of mode cutoff 
ratio have the advantage of producing relatively sim- 
ple expressions and rules of thumb. The cutoff ratio 
combines individual mode identity (eigenvalue) with 
fr«quenry and Mach number in a parameter which is 
particularly convenient for handling multimodal cases 
as a continuum. Aft fan radiation appear? to be han- 
dled well y these techniques, but the geometric 
scoustics of the inlet lip and potential flow field 
require additional work. On the other hand, the real- 
istic inlet lip geomi v snd flow field can be 
accounted for by direct numerical methods such ss the 
hybrid code. The challenge is to remove th. computer 
storage limitations on combined frequency and flow 
range. 
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The approximate analytical and the full numerical 
approaches also bracket the tange of suppressor analy- 
ses. Cutoff ratio Is again a very useful parameter 
correlating . le attenuation and optimum wall impedance 
In llt.ed ducts. The results from a full numerical 
solution, the hybrid code, have particularly empha- 
sized the Importance of modal scattering it the hard- 
aoft Interface for short treatment lengths. For the 
multimodal situation resulting from random generation 
processes, the key Input to suppressor analysis Is the 
modal energy distribution. While some broadband rases 
appear to follow equal energy. In general, the distri- 
bution must be Inferred by empirical fits to data. The 
generation models now coming Into use may provide some 
analytical guidance In this regard. The ability to 
specify acoustic treatment construction has been 
strengthened by improved understanding of the absorp- 
tion physics - * rh is concretely reflected In improved 
iipedance models particularly for Helmholtz resonator 
arrays in flow environments. The quest for wider 
absorptlor. bandwidth and enhanced low frequency atten- 
uation for a given treatment depth spurs the continued 
Investigation of bulk absorbers in full extended 
reaction configurations. 

Other Turbomachinery Installations 

Several research results highlighted in this 
review are relevant to the noise control of stationary 
turbomachinery. The Importance of reducing inflow 
disturbances to rotating blade rows, particularly for 
tone noise reduction, and the Inflow control devices 
demonstrated for turbofans apply to stationary cases. 
Tn contrast to the turbofan testing constraint, 
in-duct honeycomb devices are an option. The recent 
finding that strong tones can be generated by down- 
stream structural struts which interact with upstream 
rotors through intervening stators is another example 
of a mechanism likely to be generally operative. 

Blade pressure diagnostics represent a helpful tool 
for such mechanism Identification. Duct acoustic 
treatment can be effectively applied to pr'-du-- a 
corapcct suppressed installation. The app'oxlmate de- 
sign methods based or. cutoff ratio may be particularly 
helpful for multimodal situations that are likel to 
arise in complex, multistage machines. 
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